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ABSTRACT 

We present a spectral variability study of the XMM-Newton and Suzaku observations of one 
of the most extreme Narrow Line Seyfert 1 galaxies, IRAS 13224-3809. The X-ray spec- 
trum is characterized by two main peculiar features, i) a strong soft excess with a steep rise 
below about 1.3 keV and ii) a deep drop in flux above 8.2 keV. Although absorption-based 
interpretations may be able to explain these features by a suitable combination of ionization, 
covering factors, column densities and outflowing velocities, we focus here on a reflection- 
based interpretation which interprets both features, as well as the large soft excess, in terms of 
partially ionized reflection off the inner accretion disc. We show that the two peculiar spectral 
features mentioned above can be reproduced by two relativistic emission lines due to Fe K 
and Fe L. The lines are produced in the inner accretion disc and independently yield con- 
sistent disc parameters. We argue that the high L/K intensity ratio is broadly consistent with 
expectations from an ionized accretion disc reflection, indicating that they belong to a single 
ionized reflection component. The spectral shape. X-ray flux, and variability properties are 
very similar in the XMM-Newton and Suzaku observations, performed about 5 years apart. The 
overall X-ray spectrum and variability can be described by a simple two-component model 
comprising a steep power law continuum plus its ionised reflection off the inner accretion 
disc. In this model, a rapidly rotating Kerr black hole and a steep emissivity profile are re- 
quired to describe the data. The simultaneous detection of broad relativistic Fe L and K lines 
in IRAS 13224-3809 follows that in another extreme NLSl galaxy, IH 0707^95. Although 
the data quality for IRAS 13224-3809 does not allow us to rule out competing models as 
in IH 0707-495, we show here that our reflection-based interpretation describes in a self- 
consistent manner the available data and points towards IRAS 13224-3809 being a very close 
relative of IH 0707-495 in terms of both spectral and variability properties. These results, 
together with those based on pure broad Fe K detections, are starting to unveil the processes 
taking place in the immediate vicinity of accreting radiatively efficient black holes. 

Key words: galaxies: individual: IRAS 13224-3809 - galaxies: active - galaxies: Seyfert - 
X-rays: galaxies 



1 INTRODUCTION 

There is little doubt that narrow-line Seyfert 1 (NLSl) galaxies pro- 
vide an extreme view of the AGN phenomenon (e.g. Boiler et al. 
1996, 2002; Brandt et al. 1997; Fabian et al. 2002, 2004; Gallo 
2006; Gierlinski et al. 2008). A strong soft excess below ~2 keV 
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and significant X-ray variability are often seen, and with XMM- 
Newton came the discovery of high-energy spectral drop and cur- 
vature in some objects (e.g. Boiler et al. 2002, 2003). The physical 
interpretation is a topic of debate, and the suggestions are as ex- 
treme as the objects themselves involving relativistically blurred re- 
flection (e.g. Fabian et al. 2002, 2004; Ponti et al. 2006) or ionised 
absorption (Gierlinski & Done 2004; Middleton et al. 2007), and 
partial covering (Tanaka et al. 2004, 2005) that may be outflowing 
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(Gallo et al. 2004a). A breakthrough was the discovery of broad 
iron (Fe) La emission and reverberation-hke time lags in a long 
observation of IH 0707^95 (Fabian et al. 2009). Confirmation of 
this type of behaviour in other objects would bolster the importance 
of NLSls in understanding the innermost regions of AGN. 

IRAS 13224-3809 {z = 0.0667) along with IH 0707-495 
are perhaps the most remarkable members of this class, display- 
ing the above properties in outstanding fashion. In many ways 
IRAS 13224-3809 and IH 0707-495 are very similar objects, es- 
pecially in the X-ray band. They both show a high energy Fe K 
drop around 7-8 keV implying either an outflowing partial covering 
absorber or a large emission line (Boiler et al. 2002; Fabian et al. 
2004; Boiler et al. 2003). In both cases the material producing the 
structure must be iron overabundant with respect to the Sun. They 
also have large soft excesses (Ponti PhD thesis) with sharp drops 
around 1 keV (Leighly 1999). Moreover they show a high degree 
of variability indicating low black hole masses and high accretion 
rates. 

The motivation of this present work is to examine 
IRAS 13224-3809 in the light of the discoveries from the new ob- 
servations of IH 0707^95 (Fabian et al. 2009). To do so we make 
use of a previously analysed XMM-Newton observation (e.g. Boiler 
et al. 2003; Gallo et al. 2004b; Ponti et al. 2006) along with a new 
Suzaku observation. The purpose here is not to duplicate previous 
efforts by fitting the spectra with partial-covering (or other absorp- 
tion models) and reflection models and comparing quality of fits. 
Rather, our goal is to test the robustness of the reflection model and 
show it to be a valid interpretation, similar to that used for 1H0707- 
495. 



2 OBSERVATION AND DATA REDUCTION 
2.1 The XMM-Newton observation 

IRAS 13224-3809 was observed with XMM-Newton on 2002 Jan- 
uary 19 for about 64 ks. A detailed description of the observa- 
tion and data analysis was presented in Boiler et al. (2003) and 
Gallo et al. (2004), and we will briefly summarise details that are 
important to the analysis here. Spectral files were created from 
the original ODF using the SAS version 7.1.0. The EPIC pn and 
both MOS cameras were operated in full-frame and large-window 
mode, respectively. The MOSl and M0S2 spectra are affected by 
high background at energies greater than 6 keV, so, we did not con- 
sider MOS data in the analysis, but used these data only as a check 
of consistency. We analysed the RGS spectra, but the statistics 
hampers the detection of narrow absorption and/or emission fea- 
tures. Negligible pile up (checked with the EPATPLOT command) 
is present in the data and it does not affect the spectral results. The 
total mean EPIC-pn count rate is 2.348 cts/s, and even the brightest 
flares are well below the maximum pile-up free count rate that can 
be achived in the full-frame EPIC pn science mode (~8 cts/s; Ehle 
et al. 2003). Consistent spectral shapes are obtained whether con- 
sidering single events only or single and double events. Thus, for 
this study, we selected the latter The source plus background pho- 
tons were extracted from a circular region with a radius of 32 arc- 
sec. The background spectrum was extracted from source-free re- 
gions on the same chip as the source. With the SAS commands AR- 
FGEN and RMFGEN ancillary and response files were created. After 
filtering periods of high background (the biggest flare occurred ap- 
proximately 20 ks after the beginning of the observation and lasted 
for a few ks, see Fig.|2j we obtain net exposures of about 50 ks. The 



-1 — I — 1 — 1 — I 1 I 



— 1 1 — I — I — I — I — I — 

XMM-Newton - 



H 1 1 — h- H 



H 1 1 1 1 — h- H 



Suzciku 



_i I I I I i_ 



0.5 



—1 1— 1_J L 



Observed energy (keV) 



10 



Figure 1. Broadband .spectra of IRAS 13224-3809 observed with XMM- 
Newton pn (top panel) and Suzaku XIS (lower panel). The instrumental re- 
sponse is accounted for and the axes are identical in the two panels to ease 
direct comparison. The data are binned for appearance. 



total observed flux is 5x 10^^'^ ergs cm~^ s~^ and 2.3 x 10"^"^ ergs 
cm^^ s^^ in the 2-10 keV and 0.5-10 keV band, respectively. The 
count rate for the pn detector is 2.348 and 0.022 cts/s for the source 
and background. In our analysis the pn data has been considered in 
the 0.3-10 keV band. 

The Optical Monitor (OM) data have been analysed starting 
from the Pipeline Products (PPS). Seven and 15 exposures were 
taken through the UVW2 filter and the UV Grism, respectively. 
We only present here the analysis of the images observed with the 
UVW2 filter because of the low signal to noise of the grating ob- 
servations. The OM UVW2 flux has been corrected for Milky Way 
(Galactic) reddening. The UV data have been converted into XSPEC 
spectra files using FLX2XPS utility, part of the FTOOLS package. 



2.2 Tlie Suzaku observation 

We observed IRAS 13224-3809 with the Suzaku X-ray Telescope. 
The 200 ks run started on 2007 January 26 at 05:48:34 and ended 
on 2007 January 26 at 02:02:19 (UT). During the time the two 
front-illuminated (FI) XIS detectors (XISO and XIS3) and the sin- 
gle back-illuminated (BI) CCD (XISl) operated in standard point- 
ing mode with the target centred on the detectors (XIS nominal 
pointing). Event files were extracted from version 2. 1. 6. 15 pipeline 
processing, using a large region of 260 arcsec, as suggested by 
the ABC Guidqj, in order to prevent calibration problems. Back- 
ground counts were extracted from annuli between 345 and 470 
arcsec. Spectra and light curves were generated using XSELECT. 
Response matrices and auxiliary response files were created us- 
ing XISRMFGEN and XISSIMARFGEN (version 2007-05-14), re- 
spectively. The spectra from the two FI CCDs were examined sep- 
arately and found to be consistent with each other within uncertain- 
ties. Consequently, the XISO and XIS3 data were merged to create a 
single FI spectrum. The HXD-PIN data were processed in the stan- 
dard way, but the source was not detected by the instrument. The 



^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Figure 2. Broadband light curves of IRAS 13224-3809 observed with XMM-Newton pn in the 0.2 - 10 keV band (left panel) and with the Suzaku Fl-XlS in 
the 0.3 — 10 keV band (right panel). The time intervals used for the time resolved spectral analysis are also shown. 



total observed flux is 7x 10"^'' ergs cm~^ s~^ and 2.4x 10~^^ ergs 
cm^^ s^^ in the 2-10 keV and 0.5-10 keV band, respectively. The 
count rate for the XISO, XISl and XIS3 detectors are 0.156 cts/s 
(0.037), 0.299 cts/s (0.089) and 0.144 cts/s (0.031) for the source 
(background), respectively. In our analysis we considered the en- 
ergy band 0.4-10 keV for the XISl and 0.7-10 keV for the summed 
XISO plus XIS3. Above and below this range, deviations from rea- 
sonable models are found, indicative of calibration uncertainties. 

2.3 Joint analysis 

In the subsequent spectral analysis, source spectra were grouped 
such that each bin contained at least 20 counts and errors are quoted 
at the 90 per cent confidence level (Ax^=2.7 for one interesting pa- 
rameter). Spectral fitting was performed using XSPEC 11.3.0. All 
parameters are reported in the rest frame of the source unless spec- 
ified otherwise and all fits include Galactic absorption with the col- 
umn density fixed at its nominal value (5.3 x 10^" cm~^; Kalberla 
et al. 2005). 



3 A COMPARISON OF THE XMM-Newton AND Suzaku 
DATA 

The XMM-Newton and Suzaku observations are separated by nearly 
five years, but nevertheless appear very similar both in terms of 
brightness and shape (see Fig. [Hi. This is interesting considering 
the extreme spectral variability this object experiences on short 
(hourly) time scales (e.g. Fig.|2]l. The XMM-Newton data show a 
clear feature in the Fe K band that can be reproduced by an ab- 
sorption edge with energy E= S.IQIq'jq keV and r — l.S^J g 
(Boiler et al. 2003). The sharp drop above 8 keV is difficult to 
detect in the Suzaku spectra as the source starts to become back- 
ground dominated at those energies, however there are indications 
of it (r — LO^Q g, Ax^ ~ 12.5 for the addition of 2 new parame- 
terfl see also lower panel of Figure[Tll. Of certainty is that the fea- 

2 We observe that the drop is better detected (r = O.T^Qg, Ax^ = 16.9 
for the addition of 2 new parameters) when the spectra are extracted from a 
circular region with 2 arcmin radius (smaller than the standard 4.3 arcmin 
radius), in this way reducing the impact of the background. Nevertheless, in 
order to avoid calibration problems of the instrument effective area, in the 



ture has not shifted strongly downward in energy. The marginally 
detected 6.8 keV emission line seen in the XMM-Newton observa- 
tion (Boiler et al. 2003) is not detected in the Suzaku spectrum with 
an upper limit on its intensity of 6.7 x 10^^ ph cm^"^ s^^ and on 
the equivalent width of 1 10 eV. 

While the effective exposure was only three times longer dur- 
ing the Suzaku observation, the duration of the observations was 
substantially longer (over four days compared to three-quarters of 
a day for the XMM-Newton observation). Consequently, there was 
a possibility to examine variability on different time scales. Despite 
this, the XMM-Newton and Suzaku light curves look similar (Fig- 
ure [2j, displaying persistent variability with the largest outbursts 
being by a factor of ^ 10. In addition, the shapes of the rms spec- 
tra (see Figure|3]( are nearly identical with about a constant amount 
of variability in the 0.2-0.7 keV band, a peak between 1 and 2 keV 
and lower variability at higher energies, with a drop around 5-6 keV. 
This does suggest that the processes dominating the variability on 
short (hourly) time scales are likely the same as those which dom- 
inate on daily time scales. The main difference between the two 
rms spectra is tied to the amplitude of the drops. During the XMM- 
Newton observation the variability at low energy is about half that 
at the peak, while during the Suzaku pointing the drop is only of the 
order of 20 per cent. This difference can be explained by the pres- 
ence of a stronger constant component during the XMM-Newton 
observation (see §5.2 and Fig. 11). 

Noteworthy, is that the huge amplitude (factor of 20 — 60) 
outbursts recorded by ROSAT (Boiler et al. 1996) have never been 
seen since (Dewangan et al. 2002; Gallo et al. 2004b; this paper). 



4 UV-SOFT X-RAY EMISSION: A STANDARD DISC 
BLACK BODY COMPONENT? 

We first note that a simple power law model fails to reproduce the 
X-ray source emission because of the large soft excess. We thus 
add the contribution from a disc black body emission component 
(DISKPN, Makishima et al. 1986; Gierlinski et al. 1999). Large 
residuals are still present at ~1.3 keV and between 7 and 8 keV 

following analysis we consider the results obtained with the standard 260 
arcsec extraction radius. 
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Figure 3. RMS spectrum (see Ponti et al. 2004; Vaughan et al. 2003) from 
the XMM-Newton (top panel) and Suzaku (lower panel) observations. Both 
are calculated using ~ 5 ks bins, but the duration of the Suzaku observation 
is over five times longer. The red lines show the synthetic rms obtained from 
the best-fitting time resolved spectral variability study. The simulated rms 
have been obtained starting from the best-fitting parameters with the power 
law plus ionised disc reflection model. 



(see Fig.Q. The black body component is, however, able to repro- 
duce the soft excess emission. The best fit inner disc temperatures 
(assuming an inner disc radius of 6 Tg) are Tbb = 0.148 ± 0.004 
keV and Tbb = 0.152 ± 0.007 keV for the XMM-Newton and 
Suzaku data, respectively. The best fit luminosity of the disc black 
body component is Lo.ooi-iokcV ~ 9.7 and 6.9 x lO'*'^ erg s~^. 

Figure [5] shows the XMM-Newton (EPIC-pn) and the simulta- 
neous OM UVW^ spectra. The solid line shows the un-absorbed 
disc black body component (DISKPN) best fitting the soft X-ray 
excess emission. Once that this component is extrapolated in the 
UV band, it fails in reproducing the strong source UV emission. In 
other words, if the disc black body is responsible for the soft X~ 
ray excess, the UV flux is severely underestimated and cannot be 
accounted for by the model. 

We then fitted the UV and 2-10 keV data as described 



We also compared the OM data with the HST STIS, taken from Leighly 
& Moore (2004). The two UV data-sets agree well with each other, despite 
being acquired several years apart. We also note that the STIS spectrum 
shows a flattening at lower energies, where the strong starburst component 
becomes important (Leighly 2004). 
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Figure 4. Residuals from fitting a simple power law plus disc black body 
(DISKPN,Makishima et al. 1986; Gierlinski et al. 1999) model to the XMM- 
Newton (black) and Suzaku (both FI and BI in red). The disc black body 
emission broadly reproduce the soft excess emission, but leaves strong 
residuals at ~ 1.3 keV. 
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Figure 5. XMM-Newton EPIC pn and OM UVW2 data. The solid line shows 
the de-absorbed disc black body component that fits the soft X-ray data. The 
dotted lines show the ones fitting the UV data, assuming a black hole mass 
of 10*, 10^ and 10*' Mq, respectively from bottom to top. The dashed line 
assumes a BH mass of 5.6x10® Mq (Wang et al. 2004). A unique disc 
black body component fails in reproducing both the UV and soft X-ray 
source emission. The most probable possibility is a BH mass slightly lower 
than 10^ Mq producing the UV flux and a further emission component for 
the soft excess. 



in Vasudevan & Fabian (2007; 2009j3- The black hole mass of 
IRAS13224-3809 is not well known. The fast and high degree of 
variability suggests a small BH mass (Mbh <0.7-1x10^ Mq; 
Ponti et al. in prep.), while considerations based on the high source 
bolometric luminosity and the ultraviolet emission/absorption com- 
ponents suggest a MsH from about 6 x 10® to 10** Mq (Leighly & 
Moore 2004; Leighly 2004; Wang et al. 2004). To study the impact 
of the uncertainty associated with the poor knowledge of the BH 
mass, we consider three cases with BH masses of 10^ 10^ and 10** 
Mq. The dotted lines in Fig.|5|show the best fit model reproducing 
the UV and 2-10 keV spectrcjjfor the different BH masses. If a BH 
mass of 10^ or 10** Mq is assumed, it is not possible to reproduce 
the soft excess with a standard disc black body component because 
the disc temperature is too low. On the other hand a lower BH mass, 



The authors use a DISKPN model assuming an inner radius of 6 rg . 
^ Similar results are obtained considering the power law emission in the 
2-3 and 8-10 keV bands only. 
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in order to fit the UV data, would imply an un-physically high Ed- 
dignton accretion rate (L/L_Bdd > 100 for a 10® M©, hfLEdd ^ 1 
for a 10^ Mq). 

The dashed line in Fig. [5] shows the best fit model assum- 
ing a BH mass of 5.6 x 10** M© (Wang et al. 2004). In this case 
a super-Eddington accretion is required (rj 4), however a still 
hotter black body temperature is required to fit the soft excess. 

Modification to the disc black body emission due to atomic 
processes (in particular hydrogen and helium) are an interesting 
possibility. In fact in AGN, due to the low free-free and bound- 
free opacities in the soft X-ray regime (that allow the diffusion of 
photons from large Thomson depths below the disc surface), a high 
energy tail is generated in addition to the standard UV disk black 
body shape. However, this component, having a too low intensity 
and "temperature" (the expected maximum temperature being 50- 
60 eV), can not explain the soft excess emission of IRAS 13224- 
3809 (Ross et al. 1992). 

Alternatively, as already explored by Gierlinski & Done 
(2004), a second lower temperature Comptonization medium (with 
lower energy electrons compared to the ones producing the power 
law emission above 2 keV) may fit the soft excess spectral shape. 
IRAS 13224-3809 makes no exception. Nevertheless, we obtain (fit- 
ting the soft excess with a Comptonization component COMPST in 
Xspec) an electron temperature Tsoft ~ 0.13 keV and an opti- 
cal depth r ^ 90. These values fall within the narrow range of 
values observed by Gierlinski & Done (2004) analysing the soft 
excess of a large sample of AGN. They observed clustering of all 
the measured temperatures and optical depths to a narrow range of 
values that poses serious doubts on the second Comptonisation hy- 
pothesis for the soft excess (Gierlinski & Done 2004) and brought 
the authors to exclude this hypothesis. Thus, although the model 
can explain the soft excess in any single source, it seems to fail 
when the big picture is considered. Moreover, we stress that al- 
though soft state spectra of GBH can be modelled by two media of 
Comptonised electrons, the observed optical depth, in that case, is 
definitely smaller than r ~ 90 (as required here) resulting in a far 
smoother shape of the soft excess, than the one observed here. 

The main conclusion drawn from this exercise is that a unique 
disc black body component can not reproduce both the UV and 
soft X-ray data. The most likely situation is then that a thermal 
accretion disc component is responsible for most (if not all) of the 
UV emission, while the soft X-ray emission has a different physical 
origin. 

This is in agreement with the observation that in AGN, in gen- 
eral, when the soft excess is fitted with a disc black body compo- 
nent, the resulting temperature is observed to be too high and with 
a too small scatter, compared with expectations (only in the smaller 
mass and higher accretion rate AGN the Wien tail of the disc emis- 
sion is expected to be observed in the soft X-ray band). For this 
reason many works suggest that the real origin is tied to ionised 
absorbing/reflecting material (Gierlinski & Done 2004; Crummy 
et al. 2006). Furthermore, in observations allowing us to study the 
X-ray variability of this putative thermal component, the expected 
L oc blackbody relationship is not recovered, casting further 
doubts on this interpretation (e.g. Ponti et al. 2006). 



X-ray continuum which would be then intrinsically more intense 
and significantly steeper than observed in the 2-10 keV range. As 
shown by Boiler et al. (2003), the high energy drop at about 8 keV 
can indeed be produced by a neutral outflowing absorber partially 
covering the source. Iron overabundance is required to fit the fea- 
ture. Alternatively, a partial covering ionised absorber, similar to 
the one suggested for 1H0707-495 (Done et al. 2007), may ex- 
plain the high energy feature of IRAS 13224-3809, although the 
edge seems too sharp to be produced by partially ionized gas. 

Nevertheless, we note that to fit the source spectrum below 2 
keV two additional components are needed, an emission one (either 
a disc black body or a second power law) for the soft excess and an- 
other different absorption component to fit the ~1.3 keV structure. 
We thus fit the spectrum with a power law plus a disc black body 
emission absorbed by two ionised partial covering componento 
With one absorber roughly reproducing the high-energy drop as 
in Boiler et al. (2003), while the second absorber can reproduce 
the low energy structure only if significant blue-shifting (a ~0. 15- 
0.3 c) is allowed (fitting the structure with either O VIII or Fe M 
UTA). No good fit was found with lower outflow velocities. The 
main reason is that the model tries to fit the apparent drop below 
I keV with an L-shell iron edge, requiring high iron abundances. 
However, as Fig.[6]shows, this predicts a drop between 0.7 and 0.9 
keV due to the M-shell unresolved transition array (UTA) (Fabian 
et al. 2009; Zoghbi et al. 2010). This is a blend of numerous ab- 
sorption lines arising from the photo-excitation of ions Fel -FeXVI 
mainly produced by 2p-3d transitions (Behar et al. 2001). More- 
over the absorption-dominated model still requires a blackbody to 
model the soft excess. 

As mentioned in § 1 , the purpose of this paper is not to fit the 
complex X-ray spectrum of IRAS 13224-3809 with all possible 
models and compare the best-fitting x^- Our goal is to seek the 
simplest possible solution that can explain both the spectral and 
variability properties with a minimal set of assumptions. We sus- 
pect that combining the soft blackbody and the power law with 
further layers of partial covering absorbers with different column 
densities, ionization stages, covering factors, and with unrelated 
outflowing velocities could easily provide a fair description of the 
X-ray spectrum. However, the model would be poorly constrained 
because of its complexity and would not improve our understand- 
ing of the source. We thus consider here a different interpretation 
of the X-ray properties of IRAS 13224-3809 with the goal of de- 
scribing the data with a minimal set of physically connected model 
components. 



6 SIGNATURES OF DISC REFLECTION? 

The X-ray spectrum of IH0707-495 shows similar features at ~ 1.3 
keV and in the Fe K band. Fabian et al. (2009), Zoghbi et al. (2010) 
and the results of the detailed analysis of the RGS spectra confirm 
that the origin of the I keV feature is not due to absorption (Blustin 
& Fabian 2009). On the other hand these structures and the tim- 
ing properties of the source, are consistent with being produced by 



5 IONISED ABSORPTION 

Alternatively to the disc black body emission (or cold Comptoniza- 
tion) thermal models, the entire soft excess may be attributed to 
the effect of partially ionized gas in the line of sight, absorbing the 



° We used the same grid of models as Zoghbi et al. (2010), generated 
through the XSTAR photoionisation code (Kallman et al. 1996). The mod- 
els have abundances fixed to Solar and the illuminating flux has a steep 
power law index of 3. The fitting parameters are the absorber column den- 
sity, ionisation parameter, covering factor and Fe abundance. 
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Figure 6. The source spectrum fitted with a black body plus power law 
model absorbed by two partial covering absorber components. To attribute 
the ~1.3 keV feature to Fe L absorption requires strong Fe M absoiption at 
~ 0.8 keV that is not seen in the spectra. 

relativistic Fe L and K lines. Thus, we investigate here the consis- 
tency of the source emission with the reflection interpretation and 
its similarity with 1H0707-495. 

We thus reproduce the phenomenology of the continuum 
emission as in Fabian et al. (2009) with a power law plus disc black 
body emission. This is a crude approximation of the continuum 
emission. In particular, as shown in §4, the black body component 
is not associated with the disc emission, being in this interpretation 
associated with the reflection continuum. This study is important 
to check the similarity with 1H0707-495. We observe that, as in 
the case of 1H0707-495, a significant improvement is achieved for 
all fits when two broad disc lines with energies of ~ 1 keV and 
-E > 6 keV (depending on the continuum model) are added to the 
continuurrQ. The lines are broadened by the motion of the emitting 
material on the surface of an accretion disc and by the relativistic 
effects present in the vicinity of a black hole (in Xspec Laor pro- 
file, Laor 1991). The model parameters are the emitting energy (E), 
intensity (norm), disc inner and outer radii (ri„, Xout), disc inclina- 
tion (i) and power law disc emissivity index (q). The energies are 
broadly consistent with those expected from Fe La and Fe Kq for 
a common ionisation parameter {log(£) ~ 2 — 3.5; Kallman 1995; 
Kallman et al. 1996; 2004). It is also worth noting that the high- 
and low-energy Laor profiles are consistent with identical blurring 
parameters (see Tab.lTJ. 

Fig.[7]shows the residuals to black body plus power law con- 
tinuum fits after setting the normalisation of the Laor profiles to 
zero, clearly showing the large soft and hard residuals, both remi- 
niscent of broad and skewed emission line profiles, as expected in 
the case of relativistic disc lines. The resemblance with 1H0707- 
495 is impressive (Fabian et al. 2009). We observe that the rela- 
tivistic parameters associated with the low and high energy lines 
are consistent with each other (see Tab. 1). This suggests that the 
broadening of the low-energy feature due to the different lines of 
the Fe L complex and to other soft emission lines is likely negligi- 
ble (although it is likely that these effects are one of the the causes 
of the inferred steep disc emissivity profile). 

Given the consistency of the relativistic parameters for the two 
features, to better constrain the two broad line hypothesis, we im- 
pose the same disc inclination, emissivity index and inner disc ra- 
dius. The line energies are still consistent with Fe L and K emission 

^ A Ax^ of 385 and 79 is observed when a broad Fe L and Fe K line is 
respectively added to the source spectrum. The strong statistical variation 
indicate the importance of these features in the spectrum. 
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Figure 7. Residuals from fitting a black body plus power law continuum 
with two Laor profiles then setting the normalisation of the Laor profiles to 
zero. Black (pn), red (FI) and green (BI) models are shown. This is clearly a 
model-dependent way of showing the residulas (see Fig.|4]for comparison) 
but we believe it is useful to show it here for comparison with a similar 
figure presented in Fabian et al. (2009) for the much better statistics case of 
IH 0707^95.) 



(see lower panel of Tab. 1). Moreover, the normalization of the Fe L 
and K are in the ratio 30 ± 1 1, consistent with expectations from an 
ionised disc reflection (Ross & Fabian 2005; see also §5.2.1). The 
equivalent width of the Fe L (Fe K) line is 0.2 keV (2.8 keV) and 
0.07 keV (0.5 keV) during the XMM-Newton and Suzaku observa- 
tions, respectively. We point out that the inferred value of the Fe K 
line EW is extreme. Assuming Fe solar abundance, the Fe K line 
EW with respect to its own reflection continuum is not expected 
to be much larger than about 1 keV for spectra with F 2. A 5- 
lOxSolar Fe abundance produces a factor 2-3 increase in the line 
EW (George & Fabian 1991) which is in principle able to account 
for the observed EW if the spectrum is reflection-dominated. How- 
ever, it should be pointed out that the reflection continuum in the Fe 
K band is not a simple power law and has a skewed shape dropping 
sharply at the Fe edge. Hence a simple broad line plus power law 
fit is very likely to overestimate the line EW because a large part of 
the reflection continuum is erroneously accounted for by the broad 
line. We conclude by pointing out that the numerical results and 
parameters obtained with the phenomenological model presented 
above must be taken with caution and the model is presented here 
only as a guide for the more sophisticated and realistic reflection 
models that will be discussed in the following sections. We stress 
that, as in 1H0707-495, the two main features of the X-ray spec- 
trum, the sharp drops at ~1.1-1.2 keV and ~8.2 keV (observed 
frame), can be reproduced by two emission lines with similar ratios 
in the two objects. 

We note that, even if still too high, the inclusion in the model 
of the Fe L emission significantly lowers the temperature of the 
too hot disc component required to fit the soft excess, going from 
Tss = 0.148 - 0.152 keV to Tbb = 0.118 - 0.122 keV. Ba- 
sically, the Fe L feature well reproduces the high energy part of 
the soft excess emission. The inclusion of all the lines present in a 
disc reflection spectrum may explain the whole soft X-ray emission 
without having to invoke a thermal disc component. 
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ionised disc reflection + power law: (WABS*(KDBLUR-^REFLIONZ+POWERLAW)) 



XMM-Newron 
Suzakti 



r 

2.45±0.01 
2.38±0.02 



Abunj^g 
5±0.2 



€ 

537+16 
-25 
349±18 



q 

593+0.14 

-0.1 
4 85+0.25 

-0-2 



1 035^"J 

-0 

1.235 



inc. FkiX3_ 
52±1 3x10" 
3x10" 



Table 2. Best fit results of the XMM-Newton and Suzaku data with the parametric model composed of power law plus ionised disc relativistic reflection. We 
stress again that these errors take into account only statistical uncertainties. 



XMM-New!on 




Siizciku 






r 




r 


118 ± 4 


2.7 ±0.1 


122 ± 10 


2.5 ± 0.1 


Independent disc parameters tied between Fe L and Fe K 


FeL 


FeK 


FeL 


FeK 


Inc lOtls 
norm 93 ± 5 


6.3 ± 0.5 

6.5 ± 1.1 
^ ,+0.5 
-'•^-o.l 

°' -30 

1 a+0.8 


" -0.11 
7 + 1.5 
'■•^-0.4 

"-0.2 
37lf 


6.7 ± 0.5 

7.8 ± 2.2 

°'^-20 


Di.sc parameters tied between Fe L and Fe K 


E 0.92 + 0.02 


6.68 ± 0.3 


0.93 ± 0.04 


6.5+0.2 










Ri„ 1.34 ± 0.07 




1.38 ± 0.05 




Inc 66+ 
norm 55+^ 


1 s + 0'^ 


66 ± 2 

35;r 


1 2+°-^ 
^^■^-0.3 



Table 1. Best fit results of the XMM-Newton and Suzaku data with the 
parametric model composed of power law plus disc black body emission 
(DISKPN assuming an inner disc radius of 6 rg) and two lines broadened by 
a relativistic profile (Laor et al. 1991). The two broad lines have energies 
consistent with Fe L and Fe K. The disc black body temperature and nor- 
malisation are in units of eV. The lines energies, inner radii, inclinations and 
normalisation in keV, gravitational radii, degrees and 10~^ ph cm~^ s~^ 
in the line. Even imposing the same parameters for the relativistic profiles 
the best fit energies ai'e still consistent with emission from Fe K and Fe L 
and the line ratio with the expected one (Ross & Fabian 2005). 



6.1 Self-consistent ionised disc reflection 

An ionised reflection model (Ross & Fabian 2005) modified by the 
relativistic effects occurring on an accretion disc in the proximity of 
the black hole (e.g. Fabian et al. 2004; Crummy et al. 2006; Ponti 
et al. 2006; Petrucci et al. 2007; Larsson et al. 2008; Ponti et al. 
2009) could potentially reproduce the mean spectrum. As shown in 
§4, the disc black body emission can explain the UV emission, but 
fails in reproducing the soft excess. A prominent broad Fe L line, 
alone, can fit the high energy tail of the soft excess, but still requires 
the presence of a blackbody disc component (although a slightly 
colder one). The inclusion of the emission lines associated with the 
other elements may entirely explain the soft excess emission. For 
this reason we fit the spectrum with a simple power law model plus 
ionised reflection. 

The XMM-Newton and Suzaku spectra are fitted simultane- 
ously and in a self-consistent manner. Parameters that are not ex- 
pected to change over time scales of a few years, for example: neu- 
tral absorption (but see Gallo et al. 2007); disc inclination; and iron 
abundance, are fixed between the data. All other parameters are 
free to vary. The best fit models and the residuals from such a fit 
are shown in Table|2]and in Figure[8] 
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Figure 8. Top panel: Blurred reflection plus power law model components. 
Black (pn), red (Fl) and green (BI) models are shown. Lower panel: The 
residuals from the above fit (colour code as before). 



The model fits show a disc inclination of ~ 52 degrees; an 
inner and outer radius consistent with 1.235 and 400 Vg (Ivg = 
GM/c?) respectively, and an iron abundance ^ 5 times solar (Mor- 
rison & McCammon 1983). The disc emissivity profile is compa- 
rable at both epochs (--^ 5) as is the power law continuum photon 
index (F ~ 2.4). The ionisation parameter (defined as: ^ — l2L£; 
Tarter, Tucker & Salpeter 1969, where F is the total illuminating 
flux and n is the hydrogen number density) changes only slightly 
between the two epochs: ~ 500 erg cm s~^ during the XMM- 
Newton observation compared to ~ 340 erg cm s^^ at the later 
epoch. 
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Assuming that the reflecting disc is truncated at the innermost 
stable circular orbit (ISCO) around the black hole, the inner disc 
radius depends only on the black hole spin (see i.e. Reynolds & 
Fabian 2008). The breadth of the lines in IRAS 13224-3809 requires 
an inner radius smaller than ~2 gravitational radii (90 % confi- 
dence level), indicating a spin parameter higher than a > 0.94. 

The average 0.3 — 10 keV flux at the two epochs is also com- 
parable (~ 3 X 10~^^erg cm~^ s~^). The main difference be- 
tween the two observations is the ratio between the flux of the 
direct and reflected component, with the XMM-Newton spectrum 
being almost entirely reflection dominated (this is in agreement 
with the observed line equivalent width, see §4.1). The 10 — 50 
keV band is reflection dominated with a predicted flux of ~ 5 x 
10~^'^erg cm~^ s~^, and consistent with the lack of detection in 
the HXD. The fit is reasonably good (xl/dof = 1.13/1536) 
compared with other two-component continuum models, although 
some residuals are still present (Figure[8j. 

The advantage of the above interpretation is that the Fe La and 
Fe Kof are predicted features of the reflection model. Moreover the 
model reproduces the soft excess without the need to invoke a disc 
black body (negligible contribution is due to the extrapolation in 
the soft X-ray range of the disc black body emission fitting the UV 
data). In this case the strong UV emission would be due to the disc 
emission, while the soft X-ray by ionised reflection. Even forcing 
an additional hot disc black body component to the best fit model, 
we obtain a hot disc temperature of about 0.12 keV. This seems too 
high even for being the hard tail of the disc black body emission 
due to the low free-free and bound-free opacities present in the soft 
X-ray domain (Ross et al. 1992). 



7 SOURCE SPECTRAL VARL4BILITY 

The XMM-Newton and Suzaku broadband light curves of 
IRAS 13224-3809 (see Figure |2j exhibit large amplitude count 
rate variations, typical of this source. We explore here the time- 
resolved spectroscopy on the shortest possible time scale (set by 
requiring good quality time-resolved spectra). The XMM-Newton 
and Suzaku source light curves were sliced into 10 and 13 time in- 
tervals, respectively (see Figure |2j, and spectra were created for 
each time bin. 

7.1 Phenomenological model 

In order to gain insight on the source spectral variability and, in 
particular, to detail the variations of the Fe K and L lines, we fit 
the spectra with a phenomenological model composing of a power 
law plus a low energy disc black body continuum and two promi- 
nent broad lines representing the Fe K and Fe L emission. The disc 
black body component is required in order to approximately repro- 
duce the effect of the ionised reflection continuum. The disc black 
body temperature is fixed to the best fit values from the mean spec- 
tra (T_B_B=118 and 122 eV for XMM-Newton and Suzaku, respec- 
tively). Similarly, the parameters of the relativistic line profiles and 
their emission energies (see Tab. 1) were fixed. The model repro- 
duces well the spectra at all flux levels. 

We observe that IRAS 13224-3809 shows a correlation be- 
tween the flux and the power law spectral index, with the spectrum 
becoming steeper with increasing flux. The spectral index shows a 
large variation during both observations, going from F — 1.7 to 
about 3. This behaviour has been observed in many other sources 
(see e.g. Fig. 10 of Ponti et al. 2006) and it can be: i) either intrinsic 



due to changes in the parameters of the Comptonising phase, pro- 
ducing variations of the observed power law spectral index; ii) or 
spurious being due to the variations of the absorbing material that 
is dominating the spectral shape; iii) or due to a two component 
model comprising a steep power law varying in normalisation only 
and a harder and more constant component (McHardy, Papadakis & 
Uttley 1998; Shih, Iwasawa & Fabian 2002). The lack of a strong 
narrow component in the Fe K line (the upper limit on a neutral 
line being EW ~ 30 eV both for XMM-Newton and Suzaku and 
on a ionised being EW=120lii? eV, E=6.85j:5;;i keV for XMM- 
Newton and < 110 eV for Suzaku), excludes a strong contribution 
from a constant reflection component from distant material. Thus, 
in this latter case the constant reflection must come from the inner 
disc. 

Figure |9] shows the variations in the normalisation of Fe L 
and Fe K as a function of flux (for display purposes the inten- 
sity of Fe K has been multiplied by a factor of 30). The best 
fit values for the XMM-Newton observation are shown with filled 
blue and open light blue stars for Fe L and Fe K, respectively, 
while the Suzaku ones are shown with the filled red and open 
orange squares. The intensities of the two lines are well fitted 
by linear relations during both observation^ In fact, once the 
Fe K intensity is multiplied for a factor of 30 (see Tab. 1), we 
obtain the best fit linear relation for the Fe L and Fe K lines 
are: normFeL=axFluxo.5-iofceV with ajfMM=10.3l5^;2 x lO'^ and 
normFeK=bxFluxo.5-iofeeV with bxMM=10l4 x lO'^ ph. ergs"^ 
during the XMM-Newton observation. While the relations are: 
a=23.9±1.5x lO'^ and b=20±5x lO'^ ph. ergs"^ during the Suzaku 
one. The variations of the Fe L and K lines are thus consistent 
with being the same within each observation. This suggests that 
a unique reflection component can reproduce the two features at all 
flux states during the two observations. However, during the XMM- 
Newton observation the intensities of the lines (filled blue and open 
light blue stars) are higher than during the Suzaku observation (red 
filled and orange open squares) at the same flux level. This sug- 
gests that the reflection component (that is tracked by the emission 
lines) is correlated with the continuum within the single observa- 
tions, while for the same continuum flux the XMM-Newton obser- 
vation shows a higher amount of reflection. Nevertheless, the same 
ratio between Fe L and Fe K is always maintained. 

Figure[TO]confirms this trend. In fact although we observe that 
a fit with a constant is acceptable (x^ ~ 18.8 for 22 dof), a linear 
correlation between the intensities of the two lines provides a bet- 
ter fit than a constant (Ax^=10.2 for the same dof; ~ 8.6 for 
22 dof). This indicates that the two main spectral features at low 
(~1.2 keV) and high (~8 keV) energies are consistent with having 
a common origin. The best fit coefficient of the linear relation is: 
normFeK=axnormFeL; with a=0.028±0.006 (see Fig. 8) indicat- 
ing a ratio between Fe K and Fe L of 35±6, at all flux levels. 

7.2 Ionized disc reflection 

We fitted the spectra with a physically self-consistent model in- 
cluding a power law and an ionised disc reflection component. We 

Only for the Fe K line during the Siizuku observation a fit with a constant 
gives comparable results to the fit with a linear relation (x^ = 3.8 and 3.3, 
respectively), while during the XMM-Newton pointing a fit with a linear 
relation provides a better fit Ax^ = 6.9 (for the same dof, = 9-6 
and 2.7, respectively). The Fe L is always better fitted with a linear relation 
Ax^ = 98.2 and 105 (for the same dof, x^ = 57.3 and 7.2 during the 
Suzaku and x^ = 113.0 and 7.2, during the XMM-Newton observation). 
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Figure 9. The filled blue and open light blue stars show the Fe L and 
Fe K measurements from the XMM-Newton observation, while the filled 
red and open orange squares show the measurements from the Suzaku 
(the Fe K normalisation has been multiplied for 30). Filled symbols are 
for Fe L and open ones for Fe K. The solid and dotted hnes show the 
best fit linear relations normi^ei=axFluxo.5_iofeeV in the case of Fe L 
(a=10.3ti-2XlO^; 10^4X10'' ph. ergs"!) and Fe K (a=23.9±1.5x 10^; 
20±5xl0^ ph. ergs^^), respectively. The intensities of both broad lines 
are correlated with the continuum at each flux level. At the same flux, the 
Fe lines are significantly more intense during the XMM-Newton observa- 
tion. Nevertheless, the ratio between the Fe L and Fe K is consistent with 
being constant not only during each single observation, but also between 
the two observations taken 5 years apart. This suggests a Unk between the 
two spectral features. The Fe L/Fe K ratio is consistent with reflection from 
an ionised disc (Ross & Fabian 2005). 



assume that the inner and outer disc radius; inclination; and the el- 
emental abundances do not vary and we fix them to the best fit val- 
ues of the mean XMM-Newton and Suzaku models (consistent with 
being the same). We also fix the power law spectral index. Thus, 
the only free parameters are the normalisation of the two compo- 
nents, the disc ionisation and the emissivity index. This model re- 
produces well the spectra at all flux levels. From the best fit results 
of the time resolved spectra we simulated the expected rms spectra. 
The red line in Fig. 3 shows that this decomposition reproduces the 
source spectral variations during both observations. Despite some 
discrepancies, the general rms spectral shape is well reproduced by 
this minimal two component model. Further model complexity is 
required to accurately reproduce the rms features and it is beyond 
the scope of this paper which seeks a minimal-complexity solution. 

Fig.[TT]shows the 0.5-10 keV flux of the ionised reflection ver- 
sus the power law flux in the same band. During both observations a 
fit with a straight line is superior to a constant at a 99 % confidence 
level. The linear Pearson correlation coefficient is r = 0.75 and 
r — 0.87 for the XMM-Newton and Suzaku points, respectively. 
In both cases, the spectra appear to be reflection-dominated, most 
remarkably at low flux levels. 

The correlation between the two spectral components requires 
a positive intercept in the y-axis. Physically this indicates the pres- 
ence of an ionised disc reflection component even when the power 
law flux is null. In particular the main difference between the two 
observations is that in 2002 the amount of this residual reflection 
component was higher ('-^l.l x 10^^"^ ergs cm^^ s^^ instead of 
~6x 10~^^ ergs cm~^ s~^). The presence of this difference in the 
intensity of the residual reflection component explains also the dif- 
ference in the shapes of the rms spectra (see Fig. 3). 



o 




Figure 10. Intensity of Fe L vs. Fe K. The intensities of the two lines 
are varying in a correlated way. The variations are consistent with a con- 
stant ratio between the two features. This suggests that the two lines 
are produced by a common spectral component. The best fit relation is: 
normfe7^'=axnormj?e£ with a=0.028±0.006 indicating a ratio of 35±6, 
consistent with the values expected from disc reflection (Ross & Fabian 
2005). 



Reflection-dominated spectra have been observed in many 
AGN and in particular in NLSl accreting at a high rate (Zoghbi 
et al. 2008; Schartel et al. 2007; Fabian et al. 2004; Grupe et al. 
2008). This can happen for example when disc instabilities break 
the accretion disc in rings of dense material and the primary source, 
inside these, is thus hidden from view leaving only a reflection 
component (Fabian et al. 2002). A reflection dominated spectrum is 
also expected when the nuclear source is situated only a few grav- 
itational radii from the black hole. In this conditions light bending 
effects are strong, generating (as seen at infinity) a strong reflec- 
tion and a dim primary source (Miniutti et al. 2003; Miniutti & 
Fabian 2004). In the most extreme cases (regime I in Miniutti & 
Fabian 2004) corresponding to strongly reflection-dominated spec- 
tra (which is the relevant case here) the light bending model implies 
steep emissivity profiles and correlated variability between the re- 
flection and primary (power law) components. Although the latter 
is a prediction of the model for reflection-dominated states (as op- 
posed to less reflection-dominated cases in which the model pre- 
dicts a more contant reflection), correlated variability between the 
primary and reprocessed component is what is generally expected 
in any reprocessing model and cannot be associated unequivocally 
with the specific variability model of Miniutti & Fabian (2004). 

The upper panel of Fig.[72lshows the disc reflection emissivity 
index g as a function of the 0.5-10 keV power law flux. The emis- 
sivity index appears to be weakly anti-correlated with flux (linear 
Pearson correlation coefficient of r = —0.85), which is in line with 
the light bending model predictions in which low flux states are au- 
tomatically associated with more centrally concentrated disc illu- 
mination profiles. However, even for the most extreme light bend- 
ing effects, emissivity indices as steep as 6-7 are difficult to pro- 
duce and hence somewhat suspicious. Here we point out only one 
possible reason: the relativistic blurring model we use is based on 
the Laor (1991) code and assumes limb darkening, while numeri- 
cal simulations would in fact suggest that mild limb brightening is 
more appropriate for disc reflection (Svoboda et al. 2009). In fact, 
Svoboda et al. (2009) have shown that, by applying a limb bright- 
ning model to the best-studied broad line case of MCG-6-30-15, the 
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Figure 11. Ionised disc reflection versus power law component 0.5-10 keV 
fluxes. Both data sets give a superior flt with a linear model than with a 
constant at 99 % confidence limit (dotted and dashed line for XMM-Newton 
and Suzaku data, respectively). In particular the best fit relation is consistent 
with having the same amplitude of the variations of the two components. A 
positive intercept of the y-axis is observed extrapolating the best fit relation 
at null power law fluxes. 
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Figure 13. The disc ionization as a function of the reflection component 
flux. The observed correlation provides a sanity check for our reflection 
interpretation because both quantities should respond to the same physical 
quantity, i.e. the X-ray flux in'adiating the accretion disc. 



+ 



5x10-" 10-'- 1.5x10- 

Total 0.5-10 keV flux (ergs cm-- s-') 



8 - 



5x10-" 10-'- 1.5x10-'- 

Power Law 0.5-10 keV flux (ergs cm-- s-') 

Figure 12. Upper panel The emissivity index of the blurring profile is 
plotted vs. the 0.5-10 keV power law flux. Lower panel Disc ionisation 
parameter vs. power law flux. 



inner emissivity index can be reduced from ~5.3 to ^^3.7 for simi- 
lar inclination and black hole spin. The emissivity index may thus 
be generally overestimated, especially at relatively high inclination 
angles such as the one inferred in the present case. 

The lower panel of Fig. 12 shows the disc ionisation param- 
eter as a function of the direct (power law) flux. In this case the 
correlation is very poor (r = 0.20) meaning that the variations of 
the disc ionisation parameter do not closely follow the variations 
of the source flux (a similar behaviour is present in other sources, 
i.e. MCG-6-30-15; Ballantyne et al. 2003). We point out that, if 
light bending really plays a role in shaping the general spectral and 
variability properties of the source, the observed flux cannot be di- 
rectly connected to the flux illuminating the accretion disc (which 
is the one to which the ionization should be associated with). This 
is because, due to light bending, the disc sees a different flux than 
the observer at infinity. On the other hand, both the disc ionization 
state and the reflection intensity should be directly correlated with 
the flux irradiating the accretion disc. Hence a correlation between 
reflection intensity and ionization state is generally expected, espe- 
cially if the bulk of the reflection component comes from a limited 
region (e.g. the innermost accretion disc). Such a trend is indeed 
observed in our data (see Fig.|13|l, and represents a sanity check for 
our reflection interpretation (r = 0.65). 

IRAS 13224-3809 is not the only source where the nuclear 
emission is thought to originate from within the regions where 
strong gravity effects are not negligible (Fabian & Miniutti 2005). 
Long uninterrupted exposures of such AGN have demonstrated that 
the relation between direct emission and the disc reflection compo- 
nent is complex. The cases of MCG-6-30-15 and NGC4051 show 
that the disc reflection component is approximately constant at 
medium/high fluxes and correlated with the power law at low fluxes 
where the spectrum becomes more reflection-dominated (Ponti et 
al. 2006; Larsson et al. 2007). This is the expected behaviour when 
the intrinsic luminosity variations of the X-ray source are negligible 
and most of the flux variability is due to general relativistic effects. 
A longer observation of IRAS 13224-3809 is required to fully ad- 
dress this issue. 
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Figure 14. Frequency-dependent lags between the 0.3-1 and 1^ keV 
bands. Negative lags indicate the harder flux (dominated by the power law 
continuum) changes before the softer flux (dominated by reflection, i.e. the 
iron L line). The timing properties of IRAS 13224-3809 appear very similar 
to the ones of 1H0707-495 (Fabian et al. 2009), nevertheless the shortness 
of the observation hamper a significant detection of the lags. 



8 FREQUENCY DEPENDENT LAGS 

IRAS 13224-3809 is well known to present a high degree of com- 
plex variability (see Fig.[2]and[3}- Gallo et al. (2004) studied the 
timing properties investigating the correlations and showing a com- 
plex pattern of lags and leads between the 0.3-0.8 and 3-10 keV 
band (see Fig. 4 of Gallo et al. 2004). In combination with the 
spectral information, we performed a time-lag study (of the unin- 
terrupted XMM-Newton data) separating the energy regions dom- 
inated by the power law and ionised reflection emission. Figure 
[74l shows the frequency dependent lags between the reflection- 
dominated 0.3-1 keV and the power law-dominated 1-4 keV bands 
(see Fig. 8). Negative lags indicate that the hard (power law domi- 
nated) components leads the reflection emission. The shape of the 
frequency dependent lags closely resembles the one of 1H0707- 
495, where for the first time a reverberation lag of 30 s between the 
continuum and the energy band of the Fe L reflection line is de- 
tected (Fabian et al. 2009). In the present case a significant detec- 
tion of the lag is prevented by the short exposure (60 ks, about ten 
times lower than 1H0707-495). A long uninterrupted observation 
of IRAS 13224-3809 will clarify this issue and serve to scrutinise 
the reflection dominated interpretation proposed here. 



7.2.1 Fe L/ K energies and ratio 

The ratio between Fe L and Fe K depends on many physical param- 
eters, primarily on the power law index of the illuminating source 
and on the physical properties of the reflecting material. Kallman 
(1995) simulated the expected ratio assuming a slab of optically 
thin gas with Solar abundances photo-ionised by a X-ray power 
law source with an energy index of -1.5. In such conditions the Fe 
L/ FeK ratio is expected to be about 3 for an ionisation parameter 
log{^) ~ 2 — 3. To measure the ratio in the case of an optically 
thick accretion disc with a high Fe abundance, we computed ex- 
tensive simulations of accretion disc in the same scenario used to 
compute the ionised reflection model (Ross & Fabian 2005). In par- 
ticular the primary X-ray source illuminating the accretion disc has 
a power law spectral index of F = 2.4, the disc has a Fe abundance 
5 times higher than solar and ionisation ^ = 500 erg cm s~^. The 
ratio between the intensities of Fe L and K has been calculated con- 
sidering the Fe L region extending from 0.71 keV to 1.08 keV, and 
the Fe K region to extend from 6.31 keV to 7.07 keV (where the 
smeared Fe K-edge begins). The continuum has been fitted with a 
simple straight-line continuum in these two regions and then added 
up the line photons (how much the emergent spectrum exceeded the 
continuum fit) for each of these regions. We found a Fe L to Fe K 
photon ratio of 44.6 for this model. The reason that the L line domi- 
nates so strongly in this model is that Fe XVII dominates over most 
of the outermost Thomson depth, while the K line is suppressed 
by Auger destruction during resonance trapping. We also observe 
that, varying the ionisation parameter between 200 and 600 erg cm 
s~^, the FeL/K ratio has values of a few tens. These computations 
depend on the assumption of the straight line continuum and thus 
do not provide the precise FeL/K ratio. Nevertheless, they should 
provide reasonable values to compare with the phenomenological 
broad lines model where the continuum is also a straight line. We, 
in fact, observe that these values are roughly in agreement with the 
ones obtained from the fit with broad lines (see §6). 



9 DISCUSSION 

We have studied the X-ray spectral variability of IRAS 13224- 
3809 analyzing the XMM-Newton and Suzaku data. The source 
shows very similar behaviour during the two observations taken 
5 years apart. In particular it exhibits: i) a strong soft excess with 
a sharp feature around 1.2 keV; ii) an impressive drop in flux at 
~8.2 keV (the drop has not shifted in energy and no narrow com- 
ponent of the Fe K line is detected). 

These features may be reproduced by either absorption or ion- 
ized reflection. In the absorption interpretation the two features im- 
ply at least two distinct absorption components, while in the reflec- 
tion interpretation a single reflection component can explain all the 
major spectral features. Moreover, the reflection model naturally 
accounts also for the soft excess, while additional soft (e.g. black- 
body) components must be included in the absorption-dominated 
models we have considered. In fact when the spectrum is fitted with 
two broad lines, the energies of the lines are consistent with coming 
from Fe L and Fe K, the shapes are consistent with being produced 
by the same broadening profile and the ratio between the intensity 
broadly consistent with expectations from a dominant Fe XVII con- 
tribution. This result is very similar to what has been recently re- 
ported for a similar object 1H0707-495 (Fabian et al. 2009). In that 
case the reflection interpretation is strengthened by the detection of 
a time lag between the soft band (0.3-1 keV, dominated by the Fe L 
emission) and the power law, hampered here by the shortness of the 
XMM-Newton observation. Longer uninterrupted observations are 
required to study the timing characteristics of IRAS 13224-3809 in 
detail. 

We stress, however, that ionised disc reflection not only gives a 
reasonable fit to the mean spectrum, but also reproduces the source 
spectral variability. In particular the variability can be described 
by a variable power law with constant spectral shape and a reflec- 
tion component following the power law variations, with residual 
reflection dominating the spectrum at low fluxes. During the XMM- 
Newton observation this residual ionised disc reflection component 
is more prominent even if the source flux is similar. 

The reflection component is affected by relativistic effects. 
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strongly suggesting that it originates in the innermost regions of the 
accretion disc around a rapidly rotating Kerr black hole. The best 
fit inner radius of the accretion disc is of the order of 1.2 rg with a 
90% upper limit of ~2 rg, implying an almost maximal black hole 
spin a > 0.94. The disc illumination radial profile is consistent 
with a steep power law (q^^S-?). A weak anti-correlation between 
the illumination profile and the source direct flux is observed. This 
indicates that at lower fluxes the primary emission is likely closer 
to the back hole. During the low flux states most of the reflection 
is produced within a few gravitational radii from the black hole. In 
such conditions strong relativistic effects must be at work. 

Despite being an illustrative toy model, which most likely 
represents the zeroth-order approximation of the full relativistic 
cases, the light bending model proposed by Miniutti & Fabian 
(2004) seems to capture the main characteristics of IRAS 13224- 
3809. In particular the model predicts the existence of a regime in 
which the spectrum is strongly reflection-dominated, as we observe 
here. Within this regime, the reflection and continuum components 
are expected to correlate well, while the disc emissivity is anti- 
correlated with the direct flux (the power law flux in this case). 
Although we observe both trends, we point out that correlated vari- 
ability between the primary and reprocessed component is a gen- 
eral prediction of any reflection model. However, the light bend- 
ing model associates this behaviour only with reflection-dominated 
states, while in more standard cases, the reprocessed component is 
predicted to vary less than the primary. Future X-ray observation 
catching the source in less reflection-dominated states may be im- 
portant to test the model predictions. 

Both the strong soft excess and the extreme Fe K lines, when 
fitted with a disc ionised reflection component, require an ex- 
tremely high level of reflection, near to have a reflection domi- 
nated spectrum (see Fig. 8). We measure a ratio between the ionised 
reflection to the power law flux of 1.7 and 0.75 over the 0.001- 
1000 keV band (the entire band over which the reflection model 
is computed), while the ratio is 3.7 and 1.6 in the observed 0.5-10 
keV band for the XMM-Newton and Suzaku observations, respec- 
tively. We note that the high energy part of the reflected spectrum is 
less affected by ionisation, thus we estimated the reflection fraction 
comparing the extrapolated best fit ionised reflection to power law 
flux ratio in the 20-60 keV band (that results to be about 4.9 and 
2.1) with the one expected by neutral reflection (using a PEXRAV 
component). A reflection fraction of 7 and 4 are thus estimated in 
this way. Are such high reflection fractions consistent with the pic- 
ture described in the light bending model? As shown in Figure 15 
the observed values are rather extreme, however, they are expected 
when the primary X-ray source is less than 3-4 r^ above the cen- 
tral black hole (regime I, see Miniutti & Fabian 2004). We esti- 
mate that, in this condition, only about 10-20 % of the primary 
X-ray source is detected at infinity in the form of the power law 
emission. Thus, the measured 2-10 keV power law luminosity of 
2 — 4 X 10''^ erg s^^ would correspond to an intrinsic luminosity 
(in the frame of the source) of the order of a few times 10*"^ erg s^ ^ . 
This value is not too large if compared to the bolometric luminosity 
(Lb ~ exlO^^ergs-i), even assuming a rather extreme bolomet- 
ric correction (Elvis et al. 1995; Vasudevan et al. 2007; 2009). Also 
in agreement with such a scenario, we observe that the reflection 
fraction and the broadness of the line increase lowering the flux, 
passing from the Suzaku to the XMM-Newton observations. 

We also observe a poor correlation between disc ionization 
and observed power law flux. In the framework of the light bend- 
ing model, this is not highly surprising since the disc sees a differ- 
ent continuum than the observer precisely because of light bending. 
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Figure 15. Expected reflection fraction as a function of the source height 
above the black hole assuming a on-axis X-ray primary source including 
light bending (Fabian & Miniutti 2004). The original light bending model 
only considers continuum and Fe line production, not the reflection contin- 
uum. Thus the reflection fraction shown here has been computed by rescal- 
ing the calculated Fe line equivalent width (EW) with respect to the con- 
tinuum only using the fact that EW oc R, and by assuming that at large 
hights, where Ught bending is negligible i? = 1, as standard. 



However, the disc ionization and the reflection intensity track each 
other reasonably well, which is consistent with our reflection inter- 
pretation since they both respond to the flux irradiating the disc. 

The strong ionised reflection interpretation requires a very 
high value of the disc illumination parameter, but also of the iron 
abundance (^--^5 times solar). Nevertheless it is important to stress 
that an even higher iron abundance is required in the partial cover- 
ing interpretation (Boiler et al. 2003). Moreover supersolar metal 
abundances of the environment around supermassive black holes 
seem to be the rule, more than the exception. For example, the 
quasar environments are metal-rich with typically gas abundances 
several times Solar at all redshifts (Hamann et al. 2007). In fact, 
several methods based on optical-UV broad and narrow emission 
and absorption lines agree in observing supersolar metal abun- 
dances in large samples of quasars. The typical metallicity is ob- 
served to be roughly 2-4 times Solar (Dietrich et al. 2003a,b; Nagao 
et al. 2006; Groves et al. 2006; D'Odorica et al. 2004), in agreement 
with the value observed in IRAS 1 3224-3809. Our Galactic Centre, 
even if it is not hosting, at the moment (Ponti et al. 2010), an ac- 
creting black hole, contains large molecular clouds with Fe abun- 
dances higher than Solar. The Fe K line emitted by Sgr B2, for ex- 
ample, suggest a Fe over abundances of the order of 2, maybe more 
(Revnivtsev et al. 2004; Terrier et al. 2010). Abundances studies in 
samples of NLSl suggest that the metals are about 2-3 times higher 
than Solar in these sources (ses Nagao et al. 2002). In particular, 
a five times Solar metallicity has been measured in other nearby 
Narrow Line Seyfert galaxies (i.e. Mrk 1044, Fields et al. 2005). 
The optical/UV spectrum of IRAS 13224-3809 presents some of 
the strongest Fe emission lines of the NLSl class. The Fe over- 
abundance requires a high rate of explosion of type la supemovae 
that produce much more iron (more than an order of magnitude) 
than type II ones (Nomoto et al. 1997a,b). This may happen in the 
presence of a nuclear star cluster that formed white dwarfs that, 
through close interactions, become part of close binaries, enriching 
the environment with Iron through many SN la explosions (Shara 
& Hurley 2002; Fabian et al. 2009). 
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10 CONCLUSIONS 

IRAS 13224-3809 is a remarkable source, with one of the strongest 
soft excesses observed, and two prominent features at ~1.2 and 
~8.2 keV. 

• These features can be reproduced by 2 broad lines, in particu- 
lar: i) the best fit energies of the lines are consistent with emission 
from Fe L and Fe K; ii) the same relativistic profile can explain the 
shape of both lines; iii) the ratio between the intensities of the lines 
is broadly in agreement with expectations. 

• A single ionized disc reflection model (produced in the inner 
accretion disc) can reproduce the main features of the spectrum, 
implying a iron abundance of ~5 times solar. 

• The main spectral variations can be reproduced by a steep 
power law varying in normalisation and a reflection component 
from the inner accretion disc varying in a correlated way, but with 
a residual component, dominating at low fluxes. 

• The spectral features may, in an alternative interpretation, be 
associated with absorption components. Nevertheless in this sce- 
nario the two main features require at least two physically different 
absorbing components in highly relativistic outflow and without a 
clear physical link. A high iron abundance is also required. More- 
over, absorption-dominated models over-predict Fe M absorption 
in the soft X-rays, and they still require a likley unphysical ad- 
ditional blackbody component to describe the soft excess, casting 
further doubts on the overall interpretation. 
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